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Abstract 

Cosmological limits on neutrino masses are softened, by more than a factor 2, if 
Cold Dark Matter (CDM) and Dark Energy (DE) are coupled. In turn, a neutrino 
mass yielding up to ~ 0.20 allows coupling levels /? ~ 0.15 or more, already 
easing the coincidence problem. The coupling, in fact, displaces both P{k) and Q 
spectra in a fashion opposite to neutrino mass. Estimates are obtained through a 
Fisher-matrix technique. 
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1 Introduction 



There seem to be little doubt left: at least one neutrino mass eigenstate or, 
possibly, two of them exceed ~ 0.055 eV (direct or inverse hierarchy). This 
follows solar [T] and reactor [2] neutrino experiments, yielding Am^ 2 — 8 ^ 
10~^eV^ and, namely, atmospheric [3] and accelerator beam ^ experiments 
yielding Am^s ~ 3 x 10~^eV^. 
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Cosmology is also sensitive to neutrino mass. Since 1984, Valdarnini & Bonometto 
[53 made a detailed analysis of transfer functions in cosmologies where a part 
of Dark Matter (DM) is due to massive neutrinos, so proposing mixed DM 
models, where neutrinos play an essential role in adjusting CMB (Cosmic 
Microwave Background) anisotropies and matter fluctuation spectra to fit ob- 
servations. A large deal of work on this subject took place in the Nineties; 
mixed models were widely tested, using both the linear and the non-linear 
theory. 

Hubble diagram of SNIa [6] showed then an accelerated cosmic expansion, 
while advanced data on CMB [7] and large scale structure [S] required a spa- 
tially flat cosmology with a matter density parameter f2o,m — 0.27, so that the 
gap up to unity was to be filled by a smooth non-particle component dubbed 
Dark Energy (DE). 

All that relegated neutrinos to a secondary role in shaping cosmic data while, 
by using such advanced astrophysical data, increasingly stringent limits on 
neutrino masses could be computed (see e.g. [9]), also combining cosmo logical 
and laboratory data [TU]. Moreover, data coming from future weak lensing 
surveys seems to be powerful probes of neutrino masses [TB] . 

Standard limits on neutrino masses were recently summarized by Komatsu et 
al (2008) |lT], within the WMAP5 release, and are quoted in Table [H More 
stringent but more speculative limits are suggested in [12], who make a more 
extensive use of 2dF [13] or SDSS [H] data, and in [15], by using Lya forest 
data. 

These limits, clearly, rely on implicit assumptions concerning the dark cosmic 
sector, whose knowledge still fully relies on astrophysical data, requiring two 
components characterized by state parameters w ~ and ~ —1. But the 
assumption that no energy exchange between them occurs, tested vs. data, 
leads just to coupling limits. 

In this paper we show that spectral distortions due to CDM-DE coupling and 
to neutrino mass tend to compensate. We tentatively estimate how far we can 
go, simultaneously increasing coupling and mass, by using a Fisher Matrix 
(FM) technique. On that basis we perform a preliminary exploration of the 
parameter space. 

Table 1 

Summary of the 2-0" (95% C.L.) constraints on the sum of v masses, from WMAP 
5-year and other cosmological data sets. 





w = —1 


w 




WMAP5 


< 1.3 eV 


< 


1.5 eV 


WMAP5-hBA0+SN 


< 0.67 eV 


< 


0.80 eV 
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A large deal of work dealt with the coupling option (see, e.g. [T71I18II19||20] ). 

One of its motivations is the attempt to overcome the coincidence paradox (see 
[2 Ij ) , i.e. the fact that DE becomes relevant just at the eve of structure forma- 
tion. All that makes our epoch unique and, unless one indulges to anthropic 
views, apparently requires an explanation. However, also independently from 
this conceptual issue, our very ignorance of the physics of the dark sector re- 
quires that all reasonable options consistent with basic physics and data are 
explored. 

It is also important to outline that neutrino mass limits can be softened if DE 
with a state parameter i;; < — 1 is considered [22][23] . Unfortunately, this kind 
of state equations, yielding the so-called phantom-DE, can be justified only 
making recourse to unconventional physics. 

Here, starting from dynamical DE, we shall preliminarily discuss how a CDM- 
DE coupling can lead to a context similar to phantom DE. In this framework, 
however, no unconventional physics is involved; on the contrary, thanks to 
couphng, coincidence could be eased. 

For the sake of definiteness, in this paper we use the potential 

y(0) = (A"+Vr)exp(47r0VmJ) (1) 

admitting tracker solutions. This potential has been shown to fit WMAP data 
at least as well as ACDM [2l] and takes origin within the context of SUGRA 
theories [25] . 

The plan of the paper is as follows. In Section [2] we shall review coupled DE 
(cDE) models, comparing some aspects of its physics to phantom -DE. In 
Section [3] we show some spectra of a number of cosmological models, showing 
how CDM-DE coupling and non-vanishing neutrino mass can be selected so 
to (approximately) compensate their effects. In Section H] we debate technical 
aspects of a FM approach. In Section \5\ we give the result of such approach. 
Section [6] is then devoted to a guided exploration of the parameter space, while 
in Section [7] we present our conclusions. 



2 Coupled-DE models 

The essential feature of the scalar field (j), in order that it yields DE, is its 
self-interaction through a potential V{(f)). The simplest form of possible cou- 
pling is a linear one. It can be formally obtained by performing a conformal 
transformation of Brans-Dicke theory (see e.g. [2Z1), where gravity is modified 
by adding a (pR term [R is the Ricci scalar) to the Lagrangian. 
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Interactions with baryons are constrained by observational limits on violations 
of the equivalence principle (see, e.g. [26]) No similar constraints hold for 
CDM-DE interactions. In this case, constraints will follow from cosmological 
observations. 

In the coupled DE (cDE) scenario, as for dynamical DE, a self-interacting 
scalar field yields a cosmic component which does not cluster and has neg- 
ative pressure. As a matter of fact, its energy density and pressure read 

P = Pk + V{4>) , p = Pk- V{(t>) , (2) 

where V^cj)) is the self-interaction potential and 

Pfc = 0V2a2 . (3) 

Here dots indicate differentiation in respect to r, while the background metrics 
reads 



with dX'^ = dr^ + r\d9^ + cos^d d(f ) . (4) 



If pk ^ V , the DE state parameter approaches +1 {stiff matter) so that DE 
energy density rapidly dilutes during expansion (p oc a~^). In the opposite 
case 3> Pfc, the state parameter approaches -1 and DE allows the observed 
cosmic acceleration. In cDE models, an energy transfer occurs from CDM 
to DE, so allowing DE to have a non-negligible density since the matter- 
radiation decoupling era. However, pk is then dominant and the transferred 
energy is soon diluted. A so-called 0-matter dominated period then occurs, 
when CDM density however declines more rapidly than a~^. The increase of 
then brings to approach (the Planck mass) and V{(j)) to exceed p^. DE 
dilution then stops and DE eventually exceeds CDM density. 

Within this picture, CDM and DE stress-energy tensors {T^jj'^^\ let their 
traces read T^'^''^^'^) no longer obey separate equations; although still being 



it ought then to be 



T^''') M ^ = + CT('=)0,, (6) 
^^.^^ = - CT^'^ct^^, . (7) 



When the metric is (j4]), these equations yield 



+ 2-0 + aVi = + CaVc (8) 
a ^ 

Pc + 3-pc = - Cpccj) (9) 
a 
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Pc being CDM energy density. General covariance requires C to be a constant 
or to evolve as a function of itself. Here, instead of C, we shall use the 
dimensionless parameter 



13 = (3/167r)^/2^pC . 



(10) 



Let us then define the coupling function /(0) , through the relation 



C(0) 



rflog/ 
rf0 



(11) 



so that CDM energy density scales according to 



(12) 



Then, if we set = + pc , the eq. of motion takes the (standard) form. 



as though CDM and DE were decoupled, once the effective potential V is used. 

The CDM evolution f|T2l) . implying a density decline faster than in the absence 
of coupling, together with Eq. ([5]), implying that pc+pde has the same evolution 
as in the absence of coupling, means that pde scale dependence is different from 
what would follow from the state parameter w = Pde/Pde deducible from the 
expressions ([2]). The effective behavior, obtainable by using the potential V, 
mimics a phantom-like state equation, yielding a DE density increase with a, 
as we would find for w < —1 . 

This makes reasonable to expect that neutrino mass limits can be relaxed 
in a cDE context, as they are in the presence of phantom DE. This option, 
however, does not lead to requiring unconventional physics. On the contrary, 
if we are allowed to consider fairly high (3 values, the coincidence problem is 
also eased. 



3 Some angular and linear spectra 

The point of this paper can be appreciated through the spectra in Figures 
[T] and [H We compare a model with zero coupling and zero neutrino mass 
(00-model, hereafter) with: (i) a model with 2 massive neutrinos with mass 
nil, = 0.119 eV, yielding fl^, = 0.005 (plus 1 massless neutrino); (ii) a model 
with a CDM-DE coupling /3 = 0.049; (iii) a model with both neutrino mass 
and coupling (CM-model, hereafter). All models are spatially fiat, have adi- 
mensional Hubble parameter h = 0.71, density parameters fib = 0.04, fide = 



+ 2-0 + = , 



(13) 
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T 1 1 1 1 1 [■ 




k/hMpc"* 



Fig. 1. Transfer functions in cosmologies with/without coupling and with/without 
2 massive neutrinos. Coupling and mass are selected so to yield an approximate 
balance. The functions are multiplied by k^'^, to help the reader to distinguish 
different cases. 

0.73, spectral index = 0.96 and a cosmic opacity Topt = 0.089. DE is due 
to a SUGRA potential with A = 1.1 GeV, fitting WMAP and other data at 
least as well as ACDM. The slope a of the SUGRA potential is then fixed by 
requiring that the field density today matches Qde- 

Angular and spatial spectra are computed with an extension of the program 
CAMB [2S], allowing to treat coupled DE models also in the presence of mas- 
sive neutrinos. 

DE treatment requires that 4 first order differential equations are added to the 
basic budget. In the CAMB implementation used for this paper we account 
for the dynamical evolution of DE by using the variables 0(r) and 0(r) for 
the background DE field, as well as the variables ip{T, k) and 0(r, k) for DE 
fluctuations. 

Also CDM dynamics shall be modified when coupling is considered. The equa- 
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Fig. 2. Angular anisotropy spectra for the same models of the previous figure. Due to 
intrinsic Ci oscillations, this Figure is slightly harder to read. In the lower frame we 
also give the spectral differences between 00- and CM-models. Large I oscillations 
could be further damped by a shift by 1 or 2 units along I. The dotted lines represent 
the cosmic variance interval. 

tions, reported in Appendix A, are easily obtainable, e.g., from |29j . 

Before running CAMB we need a precursor program, to determine the value 
of a consistent with the assigned Qde and A. 

Initial conditions for background variables are easily set according to known 
tracker solutions. In the presence of coupling, the tracking regime for density 
fluctuations is complex and includes different alternatives. We however found 
that, if we set = and = at an initial time rj„, or we choose their 
expressions according to an arbitrary alternative, fluctuations accommodate 
rapidly on the right tracking and are completely independent from the initial 
choice at a time r(/c) when the k scale enter the horizon, provided that rj„ -C 

T{k). 
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In the linear theory, however, DE fluctuations matter just about the horizon 
scale and are rapidly damped afterwards. 

Both / and k ranges are selected for being those physically most significant. 
At lower Z's model discrepancies essentially vanish. In the / region shown, we 
have the sequel of maxima and minima due to primeval compression waves. 
The k range covers the scale explored by deep samples, as 2dF or SDSS, up 
to k values where non-linear effects become important. 

In the plots, spectra are multiplied by suitable powers of the abscissa / or k, so 
to reduce the ordinate range. In spite of that, in the Ci plot different spectra 
are not easy to distinguish. We then plot also the ratio ACi/Ci at constant I; 
shifts would however appear even smaller if slight shifts along the / axis (by 1 
or 2 units) were performed. 

The Figures are principally meant to show that the effects of neutrino mass 
and coupling are opposite. The coupling intensity, in fact, is selected so to 
(approximately) balance neutrino masses. 

We took, however, J2 '<^u ^ 0.67 and {3 <^ 0.075 (see PT|18|19] ): each of these 
values, by itself, is within current observational limits. Accordingly, even the 
difference between thin and thick solid-line spectra cannot be appreciated 
through current data. 

In particular, let us outline how the BAO (baryonic acoustic oscillation) struc- 
ture is faithfully reproduced when passing from the 00-models (thick solid line) 
to the CM-models (thick dashed line). 



4 Fisher matrix (data and technique) 



We then aim to test how far we can go, simultaneously increasing (3 and Vty^ 
without conflicting with data. This can be estimated by using a FM analysis 



This approach allows a rapid, semi-analytic estimate of the confldence limits 
for a speciflc experiment. It assumes a reference model as the most probable 
one, i.e. as the maximum of the likelihood distribution C{^\9) of the data 
system x given the model, described by parameters 9 = (9i). Exploiting this 
hypothesis, one can approximate £ by a multivariate Gaussian distribution, 
built using its second derivatives in respect to the parameters {6i) at the 
reference model. Nevertheless, as is known, this technique is limited by the 
actual non-Gaussian behavior of data. 
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FM is nothing but the Hessian of the log-hkehhood function: 




In the hterature, cosmological models are constrained by using a large num- 
ber of observables. To our present aims we shall directly consider the spec- 
trum of matter fluctuations P{k) and the CMB angular spectra C/^^ {XY = 
TT, TE, EE). In their recent analysis, Komatsu et al (2008) made a more 
restricted use of P{k), using only BAO's, while they used SNIa Hubble dia- 
grams, so significant also for being the first signal of DE. 

Here we chose observables directly coming from the model, in the attempt 
to leave apart observational biases, focusing just on the level of sensitivity of 
possible experiments. We consider then two different experimental contexts. 
The first one assumes that CMB spectra are measured at WMAP sensitivity 
and P{k) is measured with the sensitivity of the 2dF experiment (case W). 
The second assumes PLANCK [33] sensitivity for CMB spectra and SDSS 
sensitivity for P{k) (case P). The observational features for each mission con- 
sidered are listed in Table [2] for the case of CMB experiments and the galaxy 
surveys. 

Let us now consider first the use of CMB data only and let C^"^ be the angular 
spectra of the input model, to which we must add a white noise signal, to 
obtain 



cr = cr+N, 



XY 



tXY 



with A^, 



XY 
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xyc^x exp 



/(/ + !) 



fl2 

^FWHM 

8 In 2 



(15) 



The expressions of the Fisher matrix F}j components are then obtainable 
according to the relation 



i?u ^ nxy\nr. -iiXYX'Y' ^ r^x'Y' i^n\ 

I XY,X'Y' ^ « ^ j 



Table 2 

CMB (upper table) and galaxy surveys (lower table) specifications used in the paper. 
In the lower table, scales and volumes are in Mpc/h and (Mpc/h)^, respectively. 

Mission Imax fsky OpWHM O'T CFP 

WMAP 1000 0.8 13' 260 500 
PLANCK 2500 0.8 7.1' 42 80 

Mission kmin kmax Volume 
2dF 0.02 0.1 10^ 
SDSS 0.02 0.15 0.72x10^ 
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with 



[Cove] 



(/ + l/2)Uy 



{cr? 



{cr? 



qTEqTT 
qTEqEE 



(jTEqTT (jTEqEE 1 |-^^T£;^2 _j_ (jTT(jEEi^ 



(17) 

On the contrary, when deahng with matter power spectra, we used the follow- 
ing definition for the FM 



a, 13 * 



with 



[C0Vp]„/3 ^ 5, 



V, 



a(3 



f 



(18) 



(19) 



where Vf = {27i')^/V is the volume of the fundamental cell in k space, V is 
the volume of the survey and Vs{ka) = 4:7Tk'^6k is the volume of the shell of 
width Sk centered on ka [35|36] . In Eq. (fT9l) we left aside the contribution of 
the trispectrum, because in our analysis we considered only the linear scales, 
where the trispectrum is expected to be negligible. 



The cosmological model we consider is characterized by 9 parameters: 



Ho 

As 



opt 



Log(A/GeV) 



Physical baryon density 
Physical CDM density 
Hubble constant 
Scalar fiuctuation amplitude 
Scalar spectral index 
Reionization optical depth 

Decimal logarithm of the energy scale in SUGRA potential 
CDM — DE coupling strength 
sum of neutrino masses 



We estimate the neutrino mass density parameter, Vt^h'^, converting it from 
the total neutrino mass via 



93.5 eV 



(20) 



We compute the CMB anisotropics (temperature and polarisation) power spec- 
tra and the transfer functions, used to calculate linear matter power spectrum, 
using a modified version of CAMB. Double sided numerical derivatives were 
evaluated considering a 5% stepsize, except for A, where we adopted a 5% 
stepsize on A = Log(A/GeV). 
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10-20, 



1.06 2.12 3.18 4.24 




E.myeV 

Fig. 3. 1- and 2-o conEdence levels for a 2~-massive-neutnno model, assuming that 
the true cosmology is a SUGRA model with log(A/GeV) = 1.1, while (3 = and 
Jljy ~ 0. Thick (thin ) curves show the constraints deriving from CMB and deep sam- 
ple data (from CMB data only). Solid curves refer to the (i) case (WMAP+2dF). 
Dashed curves refer to the (ii) case (PLANCK+SDSS). In the sequel we shall exam- 
ine in detail models corresponding to the points labeled a, b, c, d, e and others. The 
location of the CM-model of Figs. [T] and [2| is indicated by an open box. The two 
locations indicated by an open circle and a cross will also be considered in detail 
below. This Figure is somehow analogous to Fig. 2 in Hannestad, 2005. 

5 Fisher matrix (results) 

In Figure [3] we then report the expected 1- and 2-cr hkehhood curves on 
the J^^u^P plane, for both cases W and P. In either case we analyse the 
constraints coming just from CMB data and those arising from the joint ex- 
ploitation of CMB and deep sample data. We performed the analysis either 
assuming 3 equal mass neutrinos, or 1 massless and 2 massive neutrinos. The 
plots shown in the Figure are obtained for the latter case, but discrepancies 
are just a minor effect. 



11 



0.25 

0.2 
0.15 

0.1 
0.05 


0.02 0.021 0.1 0.12 70 75 80 



0.25 n ^ I n FT ' ri m ' i ' ' i ' ' i ' m m i ' ' v ' i ' ' ^ ' ' i 




Log,„{A.) Log,„(A/GeV) ^opi 



Fig. 4. Correlation between (3 and the other model parameters for the W case 
(WMAP+2dF); dashed lines refer to CMB data only. 






Fig. 5. Correlation between fi and the other model parameters for the P case 
(PLANCK+SDSS); dashed lines refer to CMB data only. 

The Fisher-matrix resuhs, for the W case, substantially confirms known 1- 
and 2-0" limits on /3 [T8|19f37j . yielding /? < 0.05 and /? < 0.075, respectively, 
along the P axis {i.e. with J2^u = 0). 

On the other axis, with (3 = 0, J2 seems to be more constrained than 
what we know from current limits ( < 0.35 eV vs. J2^u < 0.8 eV 

with w 7^ —1). These discrepancies can be read as an indication of the level 
or reliability that Fisher-matrix estimate can have. In particular, they may 
be partially due to the impact of using the whole P{k) information, as well 
as to the fact that the reference cosmology is SUGRA instead of ACDM. 
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Fig. 6. 1- and 2-a limits, for the P-case, assuming that the true cosmology corre- 
sponds to the points marked by cross or open circle. The dashed lines report the 
same limits around the 00-model, as shown in Fig. O similarly, the dotted line is 
the 1-a limit around the O0~model in the W-case, as shown in Fig. [21 This Fig- 
ure shows that cosmologies, comprising CDM-DE coupling and neutrino masses, 
presently compatible with the 00 option, will be easily discriminated, at the P sen- 
sitivity level. 



However, the CMB 2-a constraint we find, J^f^f < 1-65 eV, is close to the 
95% confidence hmit < 1-5 eV obtained through a full MonteCarlo 

analysis of WMAP data only, with w ^ —1. 

The likelihood plots have the expected shape. Taken at face value they yield 
upper limits j3 < 0.22 and J^^u ^ 1-05 eV , in the case W. With the value of 
Ho used here this would correspond to fl^ ~ 0.022, more than 50% of baryon 
density. 

On the contrary, in the case P, constraints are more severe, as only CM-models 
with P < 0.07 and X^'^i/ < 0.4 eV appear consistent with the 00-model, at 
the 2-0" level. These limits are close to the maximum coupling and neutrino 
mass separately admitted in the present observational constraints. 
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In Figure H] we also show the correlations between /3 and the whole set of pa- 
rameters considered, in the W case. Correlations can be considered negligible 
for the parameters A^, Ug, A, Topt ■ The correlations with the parameters Uc, 
Uh, Ho, as expected, are stronger. Figure \5\ yields analogous results for the P 
case. 

It may also be useful to consider Figure El showing that models, including 
CDM-DE coupling and neutrino masses, compatible with the 00 option at 
the W sensitivity level, at the P sensitivity level will be well discriminated 
from it and/or also between them. 



6 Exploring the parameter space 



Of course, a FM analysis gives just a basic idea of the precision with which 
current /future data can constrain our parameter set and relax current bound 
on Qiy] on the contrary, by no means a FM analysis can tell us whether data 
favor /5 = or 7^ : here, when computing the FM from WMAP+2dF-like 
data, we choose arbitrarily a fiducial model with P = 0, flu = 0, assuming it 
to have a high likelihood on the (3-VLi, plane. 

While a FM approach allows no likelihood estimate, in this Section we wish 
to provide a few examples, exploring the parameter space that FM results 
apparently allow and, in Figures [7] and [H we exhibit the spectra for a set of 
models. 

As stated in the frame of Fig. [3, the models yielding maximum CDM-DE 
coupling (and vanishing v mass) or maximum neutrino mass (and vanishing 
coupling) have the thick line spectra. Model discrepancy is enhanced by taking 
the same amplitude Ag, instead of normalizing them to the same erg. 

The setting of models a, b, c, d, e on the J2 ^u^P plane is indicated in Fig. [3J 
They are typically within l-a boundaries. The best performance, perhaps, can 
be ascribed to models d and e. Both of them yield a present hot dark matter 
density exceeding 1 % of the critical density and 5 % of the whole DM. 

One of the scopes of this Section was testing that mildly mixed DM models, 
in the presence of a significant CDM-DE coupling, are reasonably consistent 
with data. 
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Fig. 7. Spectra for a number of cosmologies. All of them are obtained setting 
Ug = 0.96 and log^^ = —8.64, so to enhance model differences. Thick lines cor- 
respond to models presently considered in agreement with data, and yielding max- 
imum values either for ^ mi, or (3. The other lines yield models corresponding to 
the points a, h, c, d, e in Fig. [3l consistent with the OO-model at the 1-a level. 

7 Conclusions 

This paper performs a first inspection on the possibihty that high u masses 
and CDM-DE couphng yield compensating distortions of matter fluctuations 
and CMB spectra. This compensation is highly effective for small masses and 
couplings, as shown in Figs. [1] and O We then address the most significant 
question concerning the limits on coupling and u masses, when simultaneously 
considered. 

This question should be carefully addressed by using MonteCarlo techniques 
and considering all available observational constraints. Unfortunately, to do 
so, we should widen the usual parameter space, by adding 3 extra degrees 
of freedom: the coupling parameter /3, neutrino mass, and the energy scale 
A in the SUGRA model (or another equivalent parameter, in the same or in 



15 




10 100 1000 

1 

Fig. 8. Spectra of CMB anisotropies for the same cosmologies of Fig. compared 
with WMAP error amphtudes. Let us remind that all of them are obtained keeping 
the same values Ug = 0.96 and log = —8.64, so to enhance model differences. The 
relative difference of the thick line models from the 00-model appears not so wide 
as for some of the other models. Among them, however, the solid and dashed line 
models seem to perform quite well. Their performance can he improved by adjusting 
the Ho value, slightly modifying Fisher matrix outputs. 

another dynamical DE potential). 

This is among the reasons that led previous authors to perform a preliminary 
test by using a Fisher matrix technique. As is known, a FM analysis assumes 
that a given cosmology has a top likelihood and explores its neighbours. We 
can imagine, e.g. that the data favor /3 7^ in association with a significant 
u mass. A FM analysis would hardly tell us that. This calls for a MCMC 
analysis that we reserve for future work. 

Taking FM outputs at face value leads to state that models with Q^, < 0.022 
and f3 < 0.22 are observationally consistent with a. fl^y = and /3 = model. 
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Fig. 9. Density parameters for CDM, hot DM (i^'s), DE and radiation in models 
with r^jy and j3 taking the values 0.005-0.049 (short dashed), 0-0.07 (long dashed), 
0.011-0.1 (dotted), 0.02-0.21 (solid), respectively. In the last case, the DE plateau, 
extending up to the equality redshift, occurring slightly above z ~ 10^, shows a DE 
density keeping ~ 1/20 of CDM. 

Models with cDE were initially considered to overcome the coincidence prob- 
lem, in the presence of DE. To achieve this aim, one has to accept that DE is 
characterized by two different scales. For a coupled SUGRA cosmology like the 
one explored here, they are about the EW and the Planck scales. No ad-hoc 
scale is then apparently introduced, but the whole framework appears some- 
how artificial and demands for a more basic scheme, to reduce its complexity. 
An example is the double-axion model [3H], which however leads to features 
different from the ones considered here. 

Still working at a phenomenological level, in Figure [9], we however show the 
scale dependence of the density parameters, for various models with different 
Y^rUy and (3. 

In the usual case, with negligible m^,, a CDM-DE coupling compatible with 
data hardly eases the coincidence problem. Such easing is represented by the 
plateau in the Vtde curve, whose proportions are then almost insignificant. This 
does not mean that /3 7^ is not to be considered among the possible degrees of 
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freedom; e.g., in [39] it is shown that a cosmology with /3 as small as ~ 0.05, if 
inspected assuming (3 = 0, can yield wrong values for some cosmic parameters, 
including Uoc- Figure [H] however indicates that, when (3 ~ 0.2 is recovered, in 
the presence of suitably massive i/'s, a significant plateau is present and DE 
density keeps at the level ~ 1-2 % of the critical density up to 2; ~ 10^. 



Acknowledgements 



Thanks are due to R, Mainini for discussions. LPLC was supported by NASA 
grant NNX07AH59G and Planck subcontract 1290790 for this work, and would 
like to thank the Physics Department G. Occhialini, for hospitality. 



References 



[1] Q. R.Ahmad et al, Phys. Rev. Lett. 89, 011301 (2002); S. N.Ahmed et al, Phys. 
Rev. Lett. 92, 181301 (2004). 

[2] K. Eguchi et al, Phys. Rev. Lett.90, 021802h (2003); T. Araki et al, Phys. 
Rev. Lett. 94, 081801 (2005). 

[3] W. W. Allison et al, Phys. Lett.B449, 137 (1999); M. Ambrosio et al, Phys. 
Lett.B517, 59 (2001). 

[4] M. H. Ahn et al, Phys. Rev. Lett.90, 041801h (2003); D. G. Michael et al, 
Phys. Rev. Lett.97, 191801 (2006). 

[5] R. Valdarnini & S. Bonometto, A&A 146, 2 235 (1985); R. Valdarnini & S. 
Bonometto, Astrophys. J. 299, L71 (1985). 

[6] A. G. Riess, R. P. Kirshner, B. P. Schmidt, S. Jha, et al, Astrophys. J. 116, 
1009 (1998); S. Perlmutter, G. Aldering, G. Goldhaber, et al, Astrophys. J. 517, 
565 (1999); A. G. Riess, L.G. Strolger, J. Tomy, et al, Astrophys. J. 607, 665 
(2004). 

[7] P. de Bernardis, P. Ade, J. Bock, et al, Nat. 404, 955 (2000); S. Padin, J. 
Cartwright, B. Mason, et al, Astrophys. J. 549, LI (2001); J. Kovac, E. Leitch, 
c. Pryke, et al, Nat. 420, 772 (2002); P. Scott, P. Carreira, K. Cleary et al, 
Mon. Not. R. Aston. Soc.341, 1076 (2003); D. Spergel, R. Bean, Dore et al, 
Astrophys. J.Suppl. 170, 377 (2007). 

[8] M. Colless, G. Dalton, S. Maddox, et al, Mon. Not. R. Aston. Soc.329, 1039 
(2001); M. Colless, B. Peterson, C. Jackson, et al. Preprint a stro-ph/ 030658l] 
J. Loveday (the SDSS collaboration), Contemp. Phys.43, 437 (2002); M. 
Tegmark, M. Blanton, M. Strauss, et al, Astrophys. J. 606, 702 (2004); J. 
Adelman-McCarthy, M. Agueros, S. Allam, et al, Astrophys. J.Suppl. 162, 
38 (2004). 



18 



[9] A. Dolgov, Phys. Rev.370, 333 (2002); O. Elgaroy & O. Lahav, New J. Phys.7, 
61t (2005); J. Lesgourgues & S. Pastor, Phys. Rev.429, 307 (2006). 

[10] G. L. Fogli et al, Phys. Rev. D 78, 033010 (2008); G. L. Fogh et al, Phys. Rev. 
D 75, 053001 (2007). 

[11] E. Komatsu et al [WMAP Collaboration], Preprint [arXiv:0803.0547] [astro-ph] . 

[12] A. Goobar, S. Hannestad, E. Mortsell & H. Tu, JCAP 0606, 019 (2006); 
A. G. Sanchez & S. Cole, Mon. Not. R. Aston. Soc, 385, 830 (2008). 

[13] W. Percival et al., Monthly Notices of the Royal Astronomical Society 327, 
1297 (2001). 

[14] D. G. York et al, Astronomical Journal 120, 1579 (2000); C. Stoughton et 
al, Astronomical Journal 126, 485 (2002); K. Abazajian et al, Astronomical 
Journal 123, 2081 (2003). 

[15] U. Seljak et al. [SDSS Collaboration], Phys. Rev.D 71, 103515 (2005); U. Seljak, 
A. Slosar and P. McDonald, JCAP 0610, 014 (2006). 

[16] K. N. Abazajian & S. Dodelson, Phys. Rev. Lett. 91, 041301 (2003); 
S. Hannestad, H. Tu & Y. Y. Y. Wong, JCAP 0606, 025 (2006); T. D. Kitching, 
A. F. Heavens, L. Verde, P. Serra k A. Melchiorri, Phys. Rev.D 77, 103008 
(2008). 

[17] Wetterich C, A&A 301, 321 (1995); Amendola L., Phys. Rev.D62, 043511 
(2000). 

[18] L. Amendola and C. Quercellini, Phys. Rev.D 68, 023514 (2003). 

[19] Maccio' A. V., Quercellini C, Mainini R., Amendola L., Bonometto S. A., Phys. 
Rev.D69, 123516 (2004). 

[20] J. Ellis, S. Kalara, K. A. Olive & C. Wetterich, Phys. Lett.B228, 264 (1989); 
M. Gasperini, F. Piazza & G. Veneziano, Phys. Rev.D65, 023508 (2002); 
D. Comelh, M. Pietroni and A. Riotto, Phys. Lett.B 571 115 (2003); L.P. 
Chimento, A.S. Jakubi, D. Pavon & W. Zimdahl, Phys. Rev.D67, 083513 (2003); 
Z. K. Guo, N. Ohta and S. Tsujikawa, Phys. Rev.D 76 023508 (2007); E. 
Abdalla, L. R. W. Abramo, L. J. Sodre and B. Wang, Preprint arXiv:0710.1198l 
[astro-ph]; M. Manera and D. F. Mota, Mon. Not. R. Aston. Soc. 371 1373 

(2006) . 

[21] J. Khoury and A. Weltman , Phys. Rev. Lett.93 (2004) 17110; S.S. Gubser 
and J.Khoury, Phys. Rev.D 70 (2004) 104001; P. Brax et al., Phys. Rev.D 70 
(2004) 123518; D. F. Mota and D. J. Shaw, Phys. Rev. Lett. 97, 151102 (2006); 
P. Brax, C. van de Bruck, A. C. Davis, D. F. Mota and D. J. Shaw, Phys. Rev.D 
76, 124034 (2007). 

[22] S. Hannestad, Phys. Rev. Lett.95, 221301 (2005). 

[23] A. De La Macorra, A. Melchiorri, P. Serra & R. Bean, Astropart. Phys. 27, 406 

(2007) . 



19 



[24 
[25; 

[26; 
[27; 
[28; 

[29 

[3o; 

[31 

[32; 

[33 

[34 

[35 

[36; 

[37; 
[38; 



L. Colombo k M. Gervasi, JCAP 10, 001 (2006). 

P. Brax & J. Martin, Phys. Lett.B468, 40 (1999); P. Brax & J. Martin, Phys. 
Rev.DGl, 103502 (2000); P. Brax, J. Martin & A. Riazuelo, Phys. Rev.D62, 
103505 (2000). 

Damour T., Gibbons G. W. & Gundlach C., 1990, Phys. Rev., L64, 123D 
Damour T. & Gundlach C., 1991, Phys.Rev., D43, 3873 

L. Amendola, Phys. Rev.D 60, 043501 (1999); V. Pettorino and C. Baccigalupi, 
Phys. Rev.D 77, 103003 (2008) 



http://www.camb.info/ 1 



L. Amendola k D. Tocchi-Valentini, Phys. Rev.D 66, 043528 (2002) 
R. A. Fisher, J. Roy. Statist. Soc. 98, 39 (1935). 

D. S. Sivia, Data analysis: A Bayesian Tutorial, Oxford Univ. Press (Oxford, 
1996). 

M. Tegmark, A. Taylor & A. Heavens, Astrophys. J., 480, 22 (1997). 

Planck Blue Book: 

http://www.rssd.esa.int/SA/PLANCK/docs/Bluebook-ESA-SCI(2005) l 

V2.pdf 

M. C. Neyrinck & I. Szapudi, Mon. Not. R. Aston. Soc, 375, L51 (2007). 

R. Scoccimarro, M. Zaldarriaga & L. Hui, Astrophys. J. ,527, 1 (1999). 

A. J. S. Hamilton, C. D. Rimes & R. Scoccimarro, Mon. Not. R. Aston. 
Soc.,371, 1188 (2006) 

G. Olivares, F. Atrio-Barandela and D. Pavon, Phys. Rev.D 77, 063513 (2008). 

R. Mainini, S.A. Bonometto, Phys. Rev. Lett.93, 121301 (2004); R. Mainini, 
L.P.L. Colombo, S.A. Bonometto, Astrophys. J. 632, 691 (2005); R. Mainini, 
S.A. Bonometto, JCAP 0709, 017 (2007). 



[39] G. La Vacca, L. P. L. Colombo, L. Vergani & S. A. Bonometto. larXiv:0804.0"285] 
[astro-ph] . 



A Fluctuation evolution in models with dynamical and coupled 
DE 



Uncoupled DE equations - Let 

$(r, A;) = 0(r) + </?(r, A;) {Al) 
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be the dynamical DE field. In eq. (Al) the background component 0, in- 
dependent from space coordinates, and its fluctuations (p are outlined. The 
dependence on spatial coordinates is also Fourier-transformed, so to have a ip 
field dependent on the time r and the wavenumber k. 

In the absence of CDM-DE coupling, the background equation 

+ 2(d/a)(/) + aV((/)) = {A2) 

is known to hold, together with the equation 

(p + 2{a/a)ip + 4>h/2 + k'^cf) + aV"(0)(/? = {A^) 

for its fluctuations. Here h is the usual variable describing the gravitational 
field due to density fiuctuations in the synchronous gauge. 

Aside of them, the equations 

Pc + 3(d/a)pe = (^4) 

5c + kvc + h/2 = Q {Aba) 
Vc + {a/a)vc = . {A5b) 
will hold, in a generic synchronous gauge, for the CDM density pc and its 
fiuctuations 5c- Here Vc is the (gauge dependent) velocity field in CDM. 

The presence of dynamical DE clearly implies that the equation fulfilled by 
the scale factor is also modified into 



d^ = 



flcyy/a^ + {^ob + ^oc)/a^ + {Pk + V)/po,cr + Pu/Po,cr ■ {A6) 



The symbols Hg, flp'y, ^oc, ^ob have their obvious meaning; po^cr is the present 
critical energy density, while 

Pk = <PV2a^ {A7) 
is the kinetic energy density of the DE field. Finally 

exp e{mi,/xTij) 

is the energy density of z/'s with mass rriy and spin states. 



ZTT^ Jo expe(m,^/xTj,) + 1 ^ 



This equation shall not be modified by the presence of coupling, as well as the 
equations for gravitational field fluctuations, to whose source (p contributes. 
We shall omit these last equations. 

Coupled DE equations - In the presence of a constant CDM-DE coupling, 
eqs. (A2)-(A5) are modified so that their r.h.s.'s no longer vanish. More in 
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detail : 

{A2)-{AA) ^ + .... = CaVc, Pc + •••• = -Cpc^ 

(A3) (p + = Ca^pA 

{Ah) 4 + = -C0 , ^'c + = -kC^p 

All other dynamical equations keep unmodified. 
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